Strong damping of phononic heat current by magnetic excitations in SrCu2 (303)2 
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Measurements of the thermal conductivity as a function of temperature and magnetic field in the 
2D dimer spin system SrCu2(B03)2 are presented. In zero magnetic field the thermal conductivity 
along and perpendicular to the magnetic planes shows a pronounced double-peak structure as a 
function of temperature. The low-temperature maximum is drastically suppressed with increasing 
magnetic field. Our quantitative analysis reveals that the heat current is due to phonons and 
that the double-peak structure arises from pronounced resonant scattering of phonons by magnetic 
excitations. 
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During the last few years the thermal conductivity (k) 
of low dimensional spin systems has attracted consid- 
erable interest [|]-^. One reason is that in these ma- 
terials a large magnetic contribution ^'""3 to the heat 
current may be present as observed e.g. for the spin lad- 
der material Sri4_a;Ca2:Cu2404i [0. Another reason is 
that the phononic heat current k^" probes the spectrum 
of magnetic excitations as well as the spin-phonon cou- 
phng [^-||. The latter is very important in some of 
the low dimensional spin systems, e.g. in the spin-Peierls 
compound CuGeOa Both issues, the magnetic con- 

tribution to the heat current as well as the interaction 
of the phonons with magnetic excitations, are to a large 
extent unexplored and not well understood. 

A material of particular interest in this context is 
SrCu2(B03)2 (SCBO) §. The Cv?+ ions form a quasi- 
2D spin system which is by virtue of the crystal geome- 
try an experimental realization of the Shastry-Sutherland 
model In SCBO the intra-dimer and inter-dimer cou- 
phngs are of magnitude Ji ~ 72 K and J 2 — 43 K, i.e. 
J2/J1 ~ 0.6 1^. As expected for this ratio SCBO has 
a dimerized singlet ground-state separated from the ex- 
cited triplet states by a finite gap A ~ 35 K as seen in 
the magnetic susceptibility | [lO| or in inelastic neutron 
scattering . From the latter it is also known that the 
triplet excitations are almost dispersionless, i.e the group 
velocity is very small, in agreement with theoretical cal- 
culations p"^-p^. Thus, a sizeable is not expected 
for SCBO making this material a natural candidate to 
study the influence of magnetic excitations on kP^ . 

In this letter we present measurements of the thermal 
conductivity of SCBO along (kq) and perpendicular (kc) 
to the 2D spin planes in a large temperature (2.5-275 K) 
and magnetic field range (0-17 T). In zero field, both 
Ka and Kc show pronounced double-peak structures as a 
function of T. For both directions the low-T maximum 
is drastically suppressed by a magnetic field. We present 
a model based on a purely phononic thermal conductiv- 
ity and explain the double-peak structure and its field 



dependence by strong damping of the phononic heat cur- 
rent due to resonant scattering of phonons by magnetic 
excitations. With the same parameters our model con- 
siderably improves the description of the magnetic field 
dependence of the specific heat of SCBO [p^ . 

For our study two samples of rectangular-bar-shaped 
form of ^ 0.6 x 1.9 x 3 mm^ with the long direction 
along the a and c axis, respectively, were cut from larger 
single crystals of SCBO grown by the traveling solvent 
fioating zone method fl^ . The thermal conductivity was 
measured with a conventional steady state method using 
differential Chromel-Au -1-0.07% Fe thermocouples cali- 
brated in magnetic fields. Typical temperature gradients 
were of the order of 0.2 K. The absolute accuracy of our 
measurements is of order ±10% because of uncertainties 
in the sample geometry whereas the relative accuracy is 
about one order of magnitude better jiot . 

We show in Fig. |l] the thermal conductivity Ka{T) 
(= Kb in tetragonal SCBO) and Kc{T) in zero magnetic 
field. For both directions pronounced double-peak struc- 
tures are observed with the low T maxima occurring at 
« 4.5 K. For Kq the high T maximum lies at « 60 K 
and for Kc at w 30K. Above 100 K we find Kc oc T'^ 
as expected for phonon heat transport pO| whereas 
follows a T^'^-^ dependence. The magnetic field depen- 
dence of the thermal conductivity is shown in Fig. ^. The 
low T maximum is suppressed strongly by a magnetic 
field for both Kq and Kc- The magnetic field dependence 
at higher temperatures is only weak (< 2%) and close 
to the relative measurement accuracy. The minimum of 
the double-peak structure and the low T maximum sys- 
tematically shift to lower temperatures with increasing 
magnetic field. 

The behavior of k, observed here is reminiscent of that 
found in other low dimensional spin systems [Q^j^. In 
particular, in CuGeOs k, also shows an almost field- 
independent maximum at 20 K and a second one 
around 5 K, which is strongly suppressed by a magnetic 
field. This has been interpreted in terms of a magnetic 
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FIG. 1. Thermal conductivity of SrCu2(B03)2 along (ko) 
and perpendicular (kc) to the magnetic planes. The inset 
shows the experimental data (o) together with "conventional" 
phononic contributions k^'' and kJ!'' (- -) obtained by switch- 
ing off the scattering of phonons by magnetic excitations (see 
text). 

and a phononic contribution to the heat current with 
maxima at different temperatures 

For several reasons, a similar interpretation is not pos- 
sible in SCBO. Firstly, due to their flat dispersion the 
single triplet excitations are not expected to contribute 
significantly to the heat current. Secondly, k™'^9 ^ 
2D magnetic system should be strongly anisotropic, as 
is observed e.g. in Sri4_2^Ca2:Cu2404i but not in 
our measurements. The lack of a strong anisotropy ex- 
cludes also an explanation in terms of multi-triplet exci- 
tations observed in SCBO above - 5 meV 0Jl|,|l|,|||. 
These excitations move much more easily than the single 
triplet excitations , but again only within the magnetic 
planes. Finally, the relation k oc cv£ (c is the specific 
heat, V the group velocity and i the mean free path) pre- 
dicts similar T dependencies of k'""^ and the magnetic 
specific heat c'""^. But measurements of c in fields up to 
12 T (|l6| show that c™°^ gives rise to a maximum of c/T 
at temperatures close to the minima of Ka and Kc for the 
respective magnetic fields (see Fig. ||). 

From these arguments we exclude a sizable magnetic 
contribution to the heat current |Q, so that k « k^''. As 
the cause of the double peak structure and its magnetic 
field dependence we propose scattering of phonons by 
magnetic excitations. The absence of a strong anisotropy 
of Ka and Kc is then easily understood, since both, the 
anisotropy of the spin-phonon coupling and/or of the 
phonon system itself are expected to be much weaker 
than that of the 2D magnetic system. In addition, this 
scenario yields a straightforward explanation why the 
minima of and Kc occur at the maxima of c/T. The 
latter is the temperature derivative of the magnetic en- 
tropy and thus directly related to the number of magnetic 
excitations which serve as scatterers for the phonons. 

For a quantitative description we fit our data by a De- 
bye model for the phononic thermal conductivity p3[ 
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FIG. 2. Thermal conductivity Kc (a) and Ka (b) on a loga- 
rithmic temperature scale for various magnetic fields (see fig- 
ure). On increasing field the low-T maximum is suppressed. 
The lines are theoretical curves calculated for the same mag- 
netic fields via Eq. (jl|). (c) Specific heat at various fields H\\c 
(the peak height decreases with increasing field). Data points 
are from Ref. ||l^; the lines result from our theory, (d) Level 
scheme used in our model: The solid lines denote the energies 
of the ground state, 1 triplet and 2 elementary triplets, re- 
spectively. The dashed lines indicate bound states with Stot 
= or 1, which are built from two elementary triplets. B is 
the binding energy and Aq is the elementary gap. As is the 
energy difference relevant for the resonance scattering accord- 
ing to Eq. (i. 
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where v is the sound velocity, Qd the Debye tempera- 
ture, uj the phonon frequency, x = tku/ksT and r {(-u,T) 
the mean lifetime of a phonon. The total scattering rate 
T^^ is given by 
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Here T; 

processes for conventional phonon heat transport, i.e. 



, and T^j^ refer to the standard relaxation 



to boundary scattering, scattering by point defects, and 



umklapp scattering, respectively 



describes the 



scattering of phonons on magnetic excitations, which wiU 
be considered in the foUowing. 

For the double-peak structure of the phononic ther- 
mal conductivity we need a scattering mechanism that is 
most active in a certain temperature range. We therefore 
consider resonant scattering p3| , ^5| -^ . The physical pic- 
ture is the following: A phonon of energy lo is absorbed, 
thereby exciting the magnetic system (under the con- 
straint of spin conservation). Then another phonon of 
the same energy is emitted and the magnetic system is 
deexcited. If excitation (or deexcitation) starts from a 
localized state like the dispersionless triplet excitations 
of SCBO there is no directional correlation between the 
absorbed and the emitted phonon. The net effect is there- 
fore the same as that of point-like defects except for the 
frequency dependence of the matrix element Mk',k{i^), 
and the scattering rate t~^]. is of the form 



(^2_A2)2 



F{T) . (3) 



Here, As is the energy splitting of the two levels involved, 
C is the overall coupling strength between phonons and 
the two-level system and F{T) contains information on 
the state of the two- level system 1 23 2^j2^ . Notably, res- 
onant scattering starting from the magnetic ground state 
is not effective! The reason is that the ground state is 
a coherent, translational invariant state with zero mo- 
mentum so that the emitted and absorbed phonons must 
have the same momentum and energy. 

With this scenario in mind we consider resonant scat- 
tering within the magnetic excitation spectrum of SCBO 
(Fig. ||d). The lowest excited state is a single triplet with 
excitation energy Aq. It is thermally excited with proba- 
bility ci = l-[l-hmiexp(-/3Ao)]-^ Here /3 = (fc^T)-! 
and mi = 1-1-2 cosh(/3/i) contains the magnetic field de- 
pendence through h — gfisH. In the resonance pro- 
cess a second triplet is excited on one of the two adja- 
cent dimers, and it combines with the first, thermally 
excited triplet to a bound triplet state so that the to- 
tal spin is conserved. The energy of the bound triplet is 
given by 2Ao — B, where B is the binding energy [ p8[ . 
The resonance energy to be used in Eq. (|3|) is therefore 
As = (2Ao — B) — Aq = Aq — B. a complete description 
of the thermal factors of the resonance process takes ac- 
count of (i) the fact that a thermally excited triplet has 
to be present and (ii) of all resonance processes involving 
the 2 adjacent dimers. This yields 



F{T) = 2ci (A^o + A^i + (A^o - Ai)^) 



(4) 



with the probabilities Aq = 1 — ci of having a dimer in 
its singlet state and A^i = ci/3 of having a dimer in the 
triplet state with a particular component. Note that 



in the resonance process the total 5*^ component is con- 
served due to spin conservation. Therefore As does not 
depend on the magnetic field. It is only ci that is field 
dependent. 

In zero field Eqs. (^,^ yield an excellent fit (not shown) 
for Aq « 35 K and As « 20 K indicating strong binding 
effects {B w 15 K). This is in good agreement with inelas- 
tic neutron scattering and ESR experiments showing the 
lowest bound triplets at « 55 K and an elementary gap 
of Aq « 35 K (see Fig. |^) ||ri|,^. However, the decrease 
of the thermal conductivity in magnetic fields computed 
by Eqs. (||,^ is much stronger than that found in our 
experiments. 

We can account for this discrepancy by considering 
a second, magnetic field independent contribution to 
the resonance scattering. It arises from the presence of 
bound singlet excitations with energy 2Ao— 25 [p|j2^] (see 
Fig. |d). Since 2B « Aq in SCBO, the energy of these 
bound singlets is close to that of single triplet excita- 
tions yielding comparable thermal populations of these 
states. Therefore bound singlets have to be treated as 
scattering centers for resonance scattering too. In this 
case a dimer adjacent to a bound singlet is excited and 
the corresponding 3 triplets combine to Stot = 0. Since 
Stot = 0, this contribution to the resonance scattering 
does not depend on a magnetic field and lowers the sen- 
sitivity of the total resonance scattering rate to a mag- 
netic field. The contribution of the bound singlets to 
the scattering rate is treated in analogy to that of the 
single triplets. The probability that a pair of adjacent 
dimers is in a state made of two triplets is given by 
C2 = 1 — (1 + 2mi exp(— /?Ao))/^, where Z is the com- 
plete partition sum of the system of two triplets. In this 
way we achieve a theoretical description comprising sin- 
gle triplets with density per dimer rii and 2-triplet states 
with density per dimer n2. Exclusion effects - a dimer 
cannot be in a single triplet state and involved in one of 
the bound states simultaneously - are accounted for by 
ni = ci(l — 4?ii — 2n2) and 712 — 2c2(l — 2ni — 3712). 

A complication arises from a repulsive interaction v 
between neighboring triplets. The energy Aq of a sin- 
gle triplet is smaller than Ji due to virtual processes 
between neighboring dimers resulting from the dimer- 
dimer coupling. If a triplet is excited on a neighboring 
dimer these processes are (in part) blocked and Aq in- 
creases. Accounting for this on a mean-field level yields 
V = R{ni + 3n2/2) with R = (J^ - Ao)/2 In the 

calculation of k one has to replace Aq Aq + pv in all 
equations (i.e. in ci, C2 and the energies given in psf), 
where p — 2 refers to single triplets and p — 3/2 to bound 
states from two triplets. 

To illustrate the validity of our model we compute the 
specific heat for various magnetic fields. For the realistic 
parameters Aq = 36 K, i? = 18 K (from Ji = 72 K 0) 
and S = 17 K we obtain the results depicted in Fig. ac. 
At low T the agreement is very good; position and height 
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of the peak are considerably improved relative to the iso- 
lated dimer model JIgI . At higher T we presume that 
the local model developed above is too simplistic since 
the dispersion of the bound states is neglected. 

The resonant scattering, decisive for the thermal con- 
ductivity, is calculated in the same model with the same 
parameters. We obtain = Cuj'^{Ai + A2) with 

. 2ni(jVo + iVi + (iVo - jVi)^)(Ao - Bf 
^1 = TTH — ^^^5 (5a) 



((Ao - S)2 - C.2)2 



(5b) 



where TVi = (ni + 3n2/2)/3 and iVo = 1 - 3iVi. Us- 
ing this together with the usual scattering rates 
(Eq. (H)) we fit the experimental Ka and Kc by Eq. (|^). 
The sound velocity {v — 7600 m/s) along the a-direction 
and Qd = 453 K are calculated from the measured 
elastic constant cn |^^. Due to the lack of experi- 
mental data for C33 we take the same velocity for the 
fit of Kc- The relaxation rate t^^^ is then obtained 
by Tf^^^ — v/L with the characteristic sample lengths 
L — 0.75 mm (kq) and 0.5 mm (nc) of the two samples. 
The point defect scattering rate is given by r 



pt 



and the phonon umklapp scattering rate is approximated 
by r"^ = UTLO^exp{eD/uT) (24). Since the gap Aq, the 



binding energy B, and the interaction R are determined 
from the specific heat there remain only four adjustable 
parameters (P, {/, u, and C). 

As shown in Fig. |^ our model yields an almost per- 
fect fit of the measured Kc (for H — Q) over the entire T 
range. Above 10 K the fit of Ka is also almost perfect, 
whereas below 10 K a slight T shift (< 1 K) between fit 
and experimental data occurs. We use P = 7.7 ■ 10"*'^ s^ 

/K), 



(7.9 • lO-'*^ g3)^ u = 2.6- 10-31 sVK (5.7 • 10 



-31 „2 



M = 5 (10), and C = 3.35 • 10^ s^ (1.15 • 10*^ s^) for the 
calculation of Ha (kc)- These values are reasonable for 
an insulator and are comparable in magnitude to those 
found e.g. in the spin ladders 1^. Different values for Ka 
and Kc are expected in an anisotropic crystal. 

Accounting for the Zeeman splitting in ci and C2, we 
calculate Ka and k^ in magnetic fields without further 
parameter adjustment. Using g « 2.07 [0, our model 
reproduces the overall influence of the magnetic field very 
well. The high T maxima remain unchanged whereas the 
low T maxima are continuously suppressed. Again the 
agreement between calculated and experimental values is 
better for Kc than for Ka- Whereas the field influence on 
Ka is slightly overestimated by our model the experimen- 
tal values of Kc are even quantitatively reproduced up to 
the highest field. 

In the calculations, one can switch off the resonant 
scattering on magnetic excitations by setting C = 0. 
kP'' (C* — 0) obtained in this way is much larger than the 
measured thermal conductivity, as shown in the inset of 
Fig. ||. As expected for resonant scattering the damp- 



ing of the phonon heat transport is most pronounced for 
T « As. Note, however, that even at room tempera- 
ture the suppression is sizeable. This strong damping 
gives further evidence for a large spin phonon coupling 
in SCBO as has been inferred before from sound wave 
anomalies [ p9[ . 

In summary, the thermal conductivity of the low di- 
mensional quantum spin system SrCu2 (603)2 has a char- 
acteristic double-peak structure with two maxima at low 
temperatures and a pronounced magnetic field depen- 
dence. A significant magnetic contribution to the heat 
current can be excluded, since the double peak structure 
is not anisotropic and because the magnetic excitations 
are (almost) localized. Our quantitative analysis in terms 
of resonant scattering of phonons on magnetic excitations 
explains the double peak structures and their magnetic 
field dependence very well and gives evidence for strong 
spin-phonon coupling. Therefore, we call for further in- 
vestigations whether the same mechanism is at work in 
other low-dimensional spin systems. 
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